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Abstract 

The lightest neutraUno is a compehing candidate to account for cold dark matter in the 
universe in supersymmetric theories with i?-parity. In the CP-invariant theory, the 
neutralino relic density can be found in accord with recent WMAP data if neutralino 
annihilation in the early universe occurs via the s-channel A funnel. In contrast, in 
the CP-noninvariant theory two heavy neutral Higgs bosons can contribute to the 
Higgs funnel mechanism significantly due to a CP-violating complex mixing between 
two heavy states, in particular, when they are almost degenerate. With a simple 
analytic and numerical analysis, we demonstrate that the CP-violating Higgs mixing 
can modify the profile of the neutralino relic density considerably in the heavy Higgs 
funnel with the neutralino mass close to half of the heavy Higgs masses. 
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The nature of the dark matter is one of the most important questions at the interface of 
particle physics and cosmology. Recently there have been big improvements in the astro- 
physical and cosmological data, most notably due to the Wilkinson microwave anisotropy 
probe (WMAP)pP and the Sloan digital sky survey (SDSS)[?. With the data we can infer 
the following 2a range for the density of cold dark matter normalized by the critical density 

0.094 < l^cDM/i^ < 0.129, (1) 

where h ^ 0.7 is the (scaled) Hubble constant in units of 100 km/sec/Mpc. Such a precise 
determination of ^cdm^^ imposes severe constraints on any model that tries to explain it. 

In supersjnnmetric theories with i?-parity the lightest supersymmetric particle 
(LSP), which is typically the lightest neutralino x? = is stable and it serves as an 
excellent cold dark matter (CDM) candidate [HE]- However, typical mSUGRA models in 
the parameter space of minimal SUSY predict much larger values for the neutralino relic 
density than the values in the range (^. Some specific mechanisms leading to strongly 
enhanced neutralino annihilation are required to produce the observed dark matter relic 
density jH] . Such an enhancement might be due to the presence of light sleptons, enhancing 
the LSP annihilation into leptons, to an accidental degeneracy of the LSP and the lighter 
stau (or stop), leading to enhanced LSP-stau (or stop) co-annihilation, to the LSP with 
significantly mixed gaugino-higgsino components, enhancing the annihilation into gauge 
bosons, or to an accidental degeneracy ^ 2m^ with large tan/3, leading to enhanced 
annihilation through an s-channel pseudoscalar A in the CP-invariant theory. 

In particular, the enhanced LSP annihilation via a A funnel in the CP-invariant case is 
due to two reasons; (i) the S'-wave amplitude for XX — > ^ is not suppressed near threshold 
while the P-wave amplitude for xX ^ H is suppressed near threshold and (ii) the total A 
decay width^ becomes large as the A bb decay mode is greatly enhanced for large tan/3. 

The generic feature of the A funnel enhancement could, however, be greatly modified due 

to the CP-violating mixing among neutral Higgs bosons as well as due to the CP-violating 

Higgs couplings to neutralino pairs in the CP-noninvariant theory [3 |Hl El ^] • In this 

work we analyze, both analytically and numerically, the impact on the LSP relic density 

by the CP-violating Higgs mixing, loop-induced at the loop level in the CP-noninvariant 

MSSM To be specific, we consider the case when two (almost) degenerate heavy 

^The decays, A — > W^W~ and A ZZ^ arc forbidden, leading to a small A width for small tan/3 
(unless the decay ^ — + <t is open). 
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neutral Higgs bosons H and A are essentially decoupled from the lightest neutral Higgs 
boson^ and their masses are very close to twice the LSP mass. 

With the lightest neutral Higgs boson decoupled, the CP-violating mixing of the two 
nearly-degenerate heavy Higgs bosons is described by a 2 x 2 complex mass matrix, com- 
posed of a real dissipative part and an imaginary absorptive part jilOj. This mixing can be 
very large, generating frequent mutual transitions inducing large CP-odd mixing effects, 
which are quantitatively described by the complex mixing parameter X: 
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where the complex off-diagonal term A^^ of the Higgs mass matrix couples two Higgs 
states. 

The Higgs masses and widths are then shifted in a characteristic pattern by the CP- 
violating mixing pUj, of which the individual shifts can be obtained by separating real and 
imaginary parts in the relations: 
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In such a non-Hermitian mixing the ket and bra mass eigenstates have to be defined sepa- 
rately: = Cia\Ha) and {Hi\ = Cia{Ha\ {i = 2,3 and Ha = H,A); C2H = cos9,C2a = 
sin 6, C^H = — sin ^ and C^a = cos 6 in terms of the complex mixing angle 6. 

As two mass eigenstates have no definite CP parity and an enlarged mass splitting, the 
profile of the LSP relic density can considerably be modified in the heavy Higgs funnel. 
For a simple analytic and numerical illustration, we consider a specific scenario within the 
CP-violating MSSM [MSSM-CP], while a more comprehensive analysis is separately given 
in a future publication. We assume the source of CP violation to be localized entirely in 
the complex stop trilinear coupling At but all the other interactions to be CP conserving.^ 

In this situation, CP violation is transmitted through stop-loop corrections to the effec- 
tive Higgs potential, generating three CP-odd complex quartic parameters. The effective 
parameters have been calculated in Ref. [S] to two-loop accuracy and, with t/t contribu- 
tions, the parameters are determined by the parameters; the SUSY scale Ms which is taken 



^This situation is naturally realized in the MSSM in the decoupling limit with Ma > 2mz |11)[I13| . 
^This assignment is compatible with the bounds from the electric dipole moment measurements |15| . 
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to be essentially the average of two stop masses-squared, the higgsino parameter /x, the 
stop trilinear parameter At and the top Yukawa coupling ht = ^/2mt/vsm(3 defined with 
the running MS mass rht and the Higgs vacuum expectation value v ~ 246 GeV. The one- 
loop improved Born Higgs mass matrix is derived from this effective Higgs potential and 
then the matrix elements are shifted to the pole-mass parameters by including dispersive 
contributions from Higgs self-energies. 

Before evaluating the impact of the complex H/A mixing on the LSP relic density in 
the heavy-Higgs funnel, we describe an approximate procedure for estimating the relic 
density |^. The LSP number density is evolved in time according to the Boltzmann 
equation. When the temperature of the Universe is higher than the LSP mass, the number 
density is simply given by its thermal-equilibrium density. However, once the temperature 
drops below the LSP mass, the number density drops exponentially. As a result, the LSP 
annihilation rate becomes smaller than the Hubble expansion rate at a certain point when 
the LSP neutralinos fall out of equilibrium and the LSP number density in a co-moving 
volume remains constant. The present LSP relic abundance is then approximately given by 

1.07 X 109 GeV-1 

where g^, = 81 is the number of relativistic degrees of freedom and MpL = 1.22 x 10^^ GeV 
is the Planck mass. And the integral J is given by 

where (erf) is the thermally averaged LSP annihilation cross section times the relative 
velocity v of two annihilating LSPs, and Xf = m^/Tf with the freeze-out temperature 
Tf ~ m^/25 for typical weak-scale numbers. We take = 25 in the following numerical 
demonstration. 

When the heavy Higgs boson masses are large and close to twice the LSP mass, the LSP 
annihilation is dominated by heavy Higgs-boson exchanges. The LSP annihilation rate can 
then be estimated with reasonable approximation by including only the s-channel heavy 
Higgs boson exchanges. In the decoupling limit the Hxx ^"^^ ^XX couplings read 

{Xl\H\xr) = iXBrnXL)* ^ -| (iVi2 - N^, tan Ow) (sin P + cosPN,,), 
{Xl\A\xr) = {xr\A\xl)* = -^tiNi2- Nntanew){smpNi3-cospNu), (6) 
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in terms of tan/5 and the neutralino mixing matrix Ni^ {i, a = 1-4) diagonalizing the neu- 
trahno mass matrix Ai]\f as N*A4i^fN'^ = Aidiag jEl- The LSP annihilation rate multiphed 
by the relative velocity v of two LSPs can be expressed as 

^v = l E v.n^^^, (7) 

^ a,b=H,A 

where the relative velocity v is taken to be 2/3 = 2y^l — 4m^/s, and the production ampli- 
tudes Va,b and the transition decay widths Tab are defined as 

1=2,3 b=H,A 

r,, = E / d'^F D{a F)D*{b ^ F), (8) 

with the Higgs propagators Hi = l/{s — Mfj. + iMiVui)- Here, P{xX ~^ H,A) are the 
XX ^ H, A production amplitudes, determined by the couplings (0), and D{H,A F) 
the H, A ^ F decay amplitudes, for any kinematically and dynamically allowed decay 
mode F. Evaluating J{xf) in Eq. (0) with the event rate (|7j) and inserting its value into 
Eq. (jH) yields the present neutralino relic density. 

Although it is possible to calculate the masses and (transition) decay widths of the heavy 
Higgs bosons fully, we estimate them in the present work with a few approximations, which 
are reliable in the Higgs decoupling limit. In general, the light Higgs boson, the fermions 
and electroweak gauge bosons, and in supersymmetric theories, gauginos, higgsinos and 
scalar states may contribute to the loops in the complex mass matrix. In the decoupling 
limit, the couplings of the heavy Higgs bosons to gauge bosons and their superpartners are 
suppressed. Assuming all the other supersymmetric particles to be suppressed either by 
couplings or by phase space, we consider only loops by the LSP neutralino, the light Higgs 
boson and the top/bottom quark for the absorptive parts as characteristic examples; loops 
from other (s)particles could be treated in the same way of course. 

In order to demonstrate the effect of the CP-violating H/A mixing on the neutralino 
relic density in the MSSM-CP numerically, we adopt a typical set of parameters", 

M5 = 0.5 TeV, = 1.0 TeV, = 2.0 TeV; tan/5 = 5, (9) 

while varying the pseudoscalar mass M^, the SU(2) gaugino mass M2, and the phase $a 

of the trilinear term At, and taking Mi ~ O.5M2. [By reparameterization of the fields, M2 

"Analyses of electric dipole moments show that the phase of /x is quite small, unless sfermions are very 
heavy |15) : therefore its phase is set zero in our numerical demonstration. 
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is set real and positive.] For such a large fi compared to M2, the LSP is almost bino-like 
and its mass is close to Mi. 



The $^ dependence of the H/A mixing parameter X and the heavy Higgs masses and are 
displayed in Figs.IH^a) and (b), respectively, for M2,a = 0.5 TeV.** The two-state system in 
the MSSM-CP shows a very sharp resonance CP-violating mixing, purely imaginary near 
$A = O.OQvr and <I>a = O.GTvr. We note that the mass shift is indeed enhanced by more 
than an order of magnitude if the CP-violating phase rises to non-zero values, reaching 
a maximal value of the mass difference ~ 24 GeV. As a result, the two mass eigenstates 
become clearly distinguishable, incorporating significant admixtures of CP-even and CP- 
odd components mutually in the wave functions. 
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Figure 1: The dependence of (a) the real (black) and imaginary (red) parts of the 
mixing parameter X and (b) the heavy Higgs boson masses, Mh^ (black) and M^^ (red). 
M2 and Ma are set to 500 GeV. Note that 3fJe/$5mX(27r - = +3?e/ - $5mX($^) and 
the masses and widths are symmetric about (^a = tt. 



The left panel of Fig. |21 shows the allowed space of the phase $a and the normalized 

mass difference {Ma — 2m^)/2m^ for the range (JJ. Here we have set M2 to 0.5 TeV and 

**With one common phase ^a, the complex mixing parameter X obeys the relation X{2tt—^a) ~ X{^a) 
so that all CP-even quantities are symmetric when switching from ^a to 2tt — ^a- Therefore we can restrict 
the discussion to the range < ^a < tt. 
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have scanned the parameter space where 450 GeV < Ma < 550 GeV and < < tt. 
The allowed region for vr < (^a ^ Svr is simply obtained by reflecting the allowed region for 
< $A ^ with respect to $a = ^r. The green strip is for the range and the blue region 
for Vlh'^ < 0.095. In the other remaining region, we have Vlh^ > 0.129. One can clearly see 
that (i) the neutralino rehc density is indeed greatly suppressed for Ma ~ 2m^ due to the 
Higgs resonances and the detailed prediction for the relic density depends strongly on the 
value of the phase $^ as well as the mass difference between Ma and 2m^. 
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Figure 2: Left panel: The allowed phase space of the CP phase ^a (^nd the normalized 
mass difference {Ma — 2m^)/2m^ for the range (OP- The green area is for the range 
but the blue area is for the enlarged range with the lower hound ignored. Right panel: The 
allowed region of the {M2,Ma) plane for the bound Qh"^ < 0.129 in the CP-invariant case 
with ^A = (a blue strip) and CP-noninvariant case with $^ = 0.557r (two green strips). 
The values of the other relevant parameters are given in the text. 



The right panel of Fig.|2lshows the allowed regions of the (Ms, Ma) plane for nh^ < 0.129 
in the CP-invariant case (one blue strip) with $^ = and in the CP-noninvariant case 
with ^A — 0.557r (two green strips). Clearly, in order to satisfy the relic density constraint, 
the LSP mass, which is approximately O.5M2, should be close to half of the Higgs masses. 
In the CP-invariant case only the CP-odd Higgs boson A is active for the Higgs funnel 
mechanism and so only one allowed strip with its width of about 20 GeV is developed. In 
contrast, in the CP-noninvariant case with $^ = 0.557r, both of the heavy Higgs bosons 
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become active for the funnel mechanism, leading to two strips; one strip is almost identical 
to the strip in the CP-invariant case, but the other is newly developed as the i^s state, 
which is purely CP-even in the CP-invariant case, has a significant CP-odd component 
due to the CP-violating Higgs mixing. The combined width of two strips is widened due 
to the enlarged mass splitting between two mass eigenstates in the CP-noninvariant case."^''" 

To summarize. We have examined the effect of the CP-violating H/A mixing on the 
LSP annihilation cross section in the Higgs decoupling limit. By a simple analysis with a 
specific parameter set (jnj we have demonstrated that the CP-violating mixing can modify 
the profile of the LSP relic density considerably in the heavy Higgs funnel with the LSP 
mass close to half of the Higgs masses. Therefore, in order to elucidate the Higgs funnel 
mechanism through high-energy experiments on the supersymmetric particles, it is neces- 
sary to determine with good accuracy the complex mixing angle between two Higgs states 
in addition to the LSP and heavy Higgs boson masses and couplings (18) . 
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